Interfacial dynamics in transport-limited dissolution 



O 

o 

< 
en 



Martin Z. Bazant 
Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139 

(Dated: February 6, 2008) 

Various model problems of "transport-limited dissolution" in two dimensions are analyzed us- 
ing time-dependent conformal maps. For diffusion-limited dissolution (reverse Laplacian growth), 
several exact solutions are discussed for the smoothing of corrugated surfaces, including the contin- 
uous analogs of "internal diffusion-limited aggregation" and "diffusion-limited erosion". A class of 
non-Laplacian, transport-limited dissolution processes are also considered, which raise the general 
question of when and where a finite solid will disappear. In a case of dissolution by advection- 
diffusion, a tilted ellipse maintains its shape during collapse, as its center of mass drifts obliquely 
away from the background fluid flow, but other initial shapes have more complicated dynamics. 
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The analysis of interfacial dynamics using conformal 
maps ("Loewner chains") is a classical subject, which is 
finding unexpected applications in physics [l|, Ij, y] ■ For 
dynamics controlled by Laplacian fields, there is a vast 
literature on continuous models of viscous fingering |j,|5{, 
and stochastic models of diffusion-limited aggregation 
(DLA^ [3|, I3 and fractal curves in critical phenom- 
ena [l|,|^. Conformal- map dynamics has also been formu- 
lated for a class of non-Laplacian growth phenomena of 
both types I3] , driven by conformally invariant transport 
processes |3|. For growth limited by advection-diffusion 
in a potential flow, the connection between continuous 
and stochastic growth patterns has been elucidated Q, 
and the continuous dynamics has also been studied in 
cases of freezing in flowing liquids [lO|, |ll|, [ij, [iJl . 

In all of these examples, the moving interface separates 
a "solid" region, where singularities in the map reside, a 
"fluid" region, where the driving transport processes oc- 
cur and the map is univalent. (In viscous fingering, these 
are the inviscid and viscous fluid regions, respectively.) 
Most attention has been paid to problems of "transport- 
limited growth" , where the solid region grows into the 
fluid region, since the dynamics is unstable and t ypi- 
cally leads to cusp singularities in finite time [ij, |l^ |la| 
(without surface tension [j|). Here, we consider vari- 
ous time-reversed problems of "transport-limited disso- 
lution" (TLD), where the solid recedes from the fluid 
region, e.g. driven by advection-diffusion in a potential 
flow. These are stable processes, so we focus on continu- 
ous dynamics, without surface tension. 

Stochastic diffusion-limited dissolution (DLD), some- 
times called "diffusion-limited erosion" (DLE) or "anti- 
DLA" , has been simulated by allowing random walkers in 
the fluid to annihilate particles of the solid upon contact 
[iSEllilim, rather than aggregating as in DLA ^. 
Outward radial DLE on a lattice, or "internal DLA" 
(IDLA), where the random walkers start at the origin 
and cause a fluid cavity to grow in an infinite solid, has 
also been studied by mathematicians, who proved that 
the asymptotic shape is a sphere in any dimension |22l | . 

We begin our analysis by summarizing some exact so- 
lutions for continuous DLD, which could help to under- 



stand fluctuations and the long-time limits of DLE and 
IDLA. Although these solutions are known in somewhat 
different forms for Laplacian growth [ij, [iJ, |ig, it is 
instructive to summarize them prior to considering prob- 
lems of TLD by advection-diffusion, to highlight the ef- 
fects of fluid flow. 

Dissolution of surface corrugation. - Let G{w,t) be 
a conformal map from the left half plane to the fluid 
region, where steady (Laplacian) diffusion with uniform 
flux at —00 drives dissolution of the solid. The interfacial 
dynamics is given by the (dimensionless) Polubarinova- 
Gahn (PG) equation |3,la|, 



Re {G'Gt) = 1, on Rew = 0. 
An exact solution has the form, 

G{w, t)^w + h{t) + log(l + a{t)e^) 
with \a{t)\ < 1/2 (real) where 
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and C — a(0)(l — a(0)2)3/2. This solution, shown in 
Fig.n describes the decay of surface corrugations in elec- 
tropolishing [20, |2l|, |23, |2J| or the displacement of an 
inviscid fluid by an immiscible viscous fluid in a Hele- 
Shaw cell y, !l9f? ] . and closely ressembles simulations of 
"mean-field DLE" (Fig. 2 of ^). 

The linear stability of a flat interface in DLD is con- 
tained in the long-time Hmit of Eqs. ©-Q- With the 
dimensions restored {w h-> wfc, G ^-> Gfc, and t h^ kvt), 
the interface shape becomes sinusoidal. 



x{y) r^vt + Ce"*"* cos{ky) 



(5) 



for a Fourier mode of wavenumber k and mean interfacial 
velocity v. We thus recover the classical result 19, 2l|, 
l2.1 | that kv is the exponential decay rate of mode k, as 
observed in experiments J2J| . 




FIG. 1: Diffusion-limited dissolution of a corrugated sur- 
face from left to right, for /i(0) — 0, a(0) = .3, i = 
0,. 25, .50,..., 2.50 inEq. ^. 




Outward dissolution of clover-like shapes. - Next we 
consider the continuous analog of IDLA: DLD in a radial 
geometry driven by a constant diffusive flux from the ori- 
gin. This could model quasi-steady melting of an infinite 
solid around a point source of heat, or the injection of a 
viscous fluid into a Hele-Shaw celL displacing an inviscid 
fluid. The radial PG equation is |3,|5|, 



Re (wg'gt) 



1, 



1, 



(6) 



where z = g{w, t) is a univalent mapping of the interior oi 
the unit disk. A tractable case is the clover-like {N — 1)- 
fold perturbation of a circle, 



7(w, t) — ai{t)w + aNit)w 



N 



(7) 



flrst analyzed by Meyer for the (time-reversed) Hele- 
Shaw problem [lj|. For ai(0) = 1 and < c = ajv(O) < 
1/A^ (real), the solution has the implicit form. 



Na% 



1 + Nc^ + 2t and af ajv 



(8) 



Since ai{t) increases to oo from ai(0) = 1 and aN{t) 
decreases to from a7v(0) = c < 1/A^ < 1, the following 
recursion converges very quickly. 



ai 



2i-fl + iVc2(l 



-2N\ 



(9) 



and yields asymptotic approximation upon recursive sub- 
stitution. An example shown in Figure [21 illustrates the 
rapid smoothing of a four-leafed clover shape. 

From Eqs. (O-© we see that the {N - l)-fold ra- 
dial perturbation decays as a power-law, ajv oc t^^'^, 
in contrast to the exponential decay of perturbations of 
a flat interface. We conjecture that the stochastic in- 
terface in IDLA is asymptotic to the continuous DLD 
solution above, up to small logarithmic fluctuations |25l |. 
where w^ is the smallest perturbation in the initial clus- 
ter shape. In general, the continuous dynamics of DLD 



FIG. 2: Outward radial diffusion-limited dissolution driven 
by a point sink at the origin for an initial four-lobed perturba- 
tion of a circle. This exact solution is given by Eqs. (??)-(??) 
with iV = 5, c = 0.15, t = 0, 0.2, 0.4, . . . , 1.0. 



may also accurately approximate the ensemble-averaged 
stochastic dynamics of IDLA, which is not the case for 
DLA and other unstable aggregation processes J9j . 

An (apparently new) explicit solution is possible for 
N = 2. In that case, Equation (|2J) reduces to a depressed 



cubic |2^, 



3pa2 — 2q, solved by the formula of Car- 



dano and dal Ferro, 



ai = 



V 9 + Vi^ - p^ + yq- V 



q2 — p3 



(10) 



where p == -(t + c^ + l/2)/4 and q = -c/4. 

Inward dissolution of star-like shapes. - Next we 
briefly consider inward DLD with a constant diffusive flux 
at inflnity, which simply is the time-reverse of Laplacian 
growth in a radial geometry Q . Now the map in Eq. ® 
must be univalent outside the unit disk, with a Laurent 
expansion. 



liw,t)= J2 a«(Ow^", \w\ >1 



(11) 



There are well-known solutions for (A^ -|- l)-fold pertur- 
bations [ilEl, 

g{w,t)=ai{t)w + a-Nit)w^^ for \w\ > 1, (12) 

where, without loss of generality, ai(0) = 1 and 
a_Ar(0) = c < 1/A^ is real. Again, the Laurent coeficients 
satisfy a pair of nonlinear equations, which is most easily 
solved as a flxed-point iteration. 



ai ^ ^1 - 2t + Nc^iaj^ - 1) 



and a_ 



N 



ca 



N 



(13) 




FIG. 3: Inward radial difFusion-limited dissolution driven 
by a sink at cx3 for a four-pointed shape. This exact so- 
lution is given by Eqs. ^-^ with N = S, c ^ 0.3, 
t = 0, 0.04, 0.08, . . . , 0.36, 0.3649; the collapse occurs at fc = 
0.365. 



The only qualitative difference with outward DLD is that 
the solid collapses to a point in a finite time, tc = (1 — 
N(?)I2. For A^ = 1, an ellipse (0 < c < 1) or circle 
(c — 0) maintains its shape during collapse, 



g{w,t) 



1 



(w+-) 



(14) 



For A^ > 1, the shape approaches a circle prior to col- 
lapse, according to the asymptotic formula. 



ai 



(i) ^ ^2{ta~t)+Nc^2it,~t))^ (15) 

with a_jv(t) ~ cai{t)^ . The collapse of a four-pointed 
shape {N — 3) is shown in Figure O 

Advection- diffusion-limited dissolution. - The dynam- 
ics of dissolution become more interesting when driven by 
non-Laplacian (but conformally invariant |3]) transport 
processes J7|, where right-hand side of the PG equation 
© is replaced by the nonuniform, time-dependent flux 
to the interface, a{w,t). An important example is dis- 
solution by advection-diffusion in a potential flow, e.g. 
the erosion of rock by flowing water, the evaporation of 
a fiber coating in a flowing gas, or the melting of a solid 
column in a flowing liquid. The time-reversed growth 
problem has been studied extensively in the contexts of 
freezing [ifl Jl] ljl2. Il3l| and advection-diffusion-limited 
aggregation 0, , but it seems that dissolution - which 
leads to collapse in finite time - has not been analyzed. 
For a given initial shape and background fiow, when and 
where will collapse occur? 

Consider a finite solid of constant concentration and 
arbitrary shape in a two dimensional potential flow of 



zero concentration and uniform velocity far away. The 
relative importance of advection to diffusion is measured 
by the Peclet number. Peg — UL/D, for a background 
fluid velocity U, diffusivity D and length L. The time- 
dependent Peclet number, Pe{t) — Peoai(i), is defined 
by the conformal radius, ai{t), in Eq. ((TT|l . When sohd 
dissolves (ai(i) — *• 0), diffusion eventually dominates, 
and it is natural to focus on the low-Pe limit. 

The flux proflle a{d, Pe) on the absorber has been stud- 
ied extensively, and very accurate asymptotic approxima- 
tions are available [23. (A numerical code in matlab is 



also at http://advection-diffusion.net,) From the low-Pe 
approximation. 



lojPe) 
KoiPe/2)' 



-Pe cos( 



Pe 



dte*^ 



,h{t) 



t 



(16) 
which is uniformly accurate in angle 9 up to Pe =10 ^, 
let us keep only the leading terms. 



l + Pe cos 9 
log(Pe/4) - 7 



Pe cos 9 



(17) 



where 7 = 0.577215 ... is Euler's constant. 

In the flnal stage of collapse where — logPe ^ 7, 
the interface is asymptotically circular, g{w,t) ^ ai{t)'w. 
From (|17|l . the radius in this regime satisfles, aidai/dt ^ 
— l/logfli, and thus has the asymptotic form. 



ai 



4(t - t,) 



log t^- log log 7^ 



(18) 



To describe the dynamics starting at small Pe and end- 
ing just prior to collapse, it is reasonable to further set 
logPe « constant. In this intermediate regime, the in- 
terfacial dynamics is given by 



Ke{wg'gt) = —l + bai[t)cos9 for 



(19) 



(with a suitable choice of units). The cos 6* term has 
the effect of exciting new modes in the shape of the dis- 
solving solid, which are not present in the initial condi- 
tion. We will see that a circle, g{w,0) — w, translates 
away from the flow, g{w, t) = w + ao(i), as its upstream 
side dissolves more quickly prior to collapse. It can also 
be shown that an {N + 1)— fold perturbation of a circle 
(|12|l will translate and develop an A^— fold perturbation, 
g{w,t) = ai{t)w + ao{t) + a_Ar+i/w^~^ -|- aAr(t)/w^. 
Higher Fourier modes in a{w,t) excite additional terms 
in the Laurent expansion of the shape, g{w,t). 

Collapse of a tilted ellipse in a uniform flow. - To 
illustrate these principles, let us consider a solid ellipse, 
g{w, 0) = w + c/w (c < 1), which is tilted at an angle (fi — 
(argc)/2 with respect to a background flow in the positive 
X direction. During dissolution, the shape remains an 
ellipse, but its center of mass, ao(i), move. An example 
is shown in Fig. ^ 




FIG. 4: Inward advection-difFusion-limited dissolution of a 
tilted elliptical solid in a uniform background flow from left 
to right (c = 0.5i, b = 0.8, t = 0, 0.04, 0.08, . . . , 0.36, 0.3749). 
The major axis of the ellipse remains oriented at 7r/4, while 
its center of mass moves at an angle of 7r/6 relative to the 
background flow. At time ic ~ 0.375, the solid disappears at 
the point ao(tc) = 0.4 + 0.2i. 



from the flow, but also away from the end of the elhpse 
which protrudes upstream. The velocity of the center of 
mass is constant and in the 1 + c direction, 



ao{t) 



b{l + c)t 
1-lcP ■ 



(21) 



Note that the center of mass does not move along the 
major axis of the ellipse at angle 0, but instead at an 
oblique angle 9 given by sinO = |c|sin(2(/) — 9). The 
final collapse occurs at the point, ao(tc) = b{l + c)/2. 
In the simplest case, c — 0, a circle maintains its shape 
while its center of mass translates away from the flow at 
(dimensionless) velocity b until collapse occurs at a; = b/2 
at time tc = b/2. 

For any initial shape, the center of mass will drift away 
from the flow, as well as away from any protrusions in 
the direction of the flow. For non-elliptical shapes, it 
is nontrivial to predict the exact time and place of col- 
lapse, even for the simplifled dynamics of Eq. p9|l . For 
more general transport-limited dynamics 7], predicting 
the collapse seems like an interesting open problem. 

The author thanks J. Propp for an introduction to 
IDLA. 



The Laurent expansion is given by g{w, t) — ai{t)w + 
ao{t) + a-i{t)/w, where ai{t) is real, but ao{t) and 
a_i(i), are complex. Substituting into H19|) and inte- 
grating shows that the conformal has a square-root sin- 
gularity. 



aiW-A/l 



t 



tc^ 



1 



\c\' 



(20) 



since the area decreases linearly to zero (in the approxi- 
mation of constant total flux): A{t) = A{0) — 27ri, where 
A{0) = 7r(l — |cp). The collapse time depends on the 
initial shape through |c|. The slowest collapse, t^. = 1/2, 
occurs for a circle, c = 0, while the collapse time tends 
to zero for a very elongated elhpse, \c\ — > 1, regardless of 
orientation. 

For a constant total flux, the collapse time tc does 
not depend on the bias introduced by the flow veloc- 
ity (through b), although the flow affects the time scale 
through the initial Peclet number, Peg. Mathematically, 
this is a general consequence of the conformal invariance 
of the transport process [3 , which causes the total flux 
(Nusselt number) to depend only on the conformal radius 
and not the asymmetric shape of the particle j27| . Phys- 
ically, the enhancement of dissolution on the upstream 
side of the solid is cancelled by the reduction in dissolu- 
tion on the downstream side. 

During dissolution, the ellipse keeps its shape and 
its orientation with respect to the flow direction since 
a_i(i) = coi(i). However, the center of mass moves away 
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